Volume contraction at the Jahn- Teller transition of LaMnOs 



o 
o 

(N 

> 

o 
O 

ON 



i 

-4-* 



i 

c 

o 
o 



> ■ 
o ■ 
o : 
r» . 
o ■ 

co : 
o . 



c : 
o 

o ■ 



X 



T. Maitra, 1 P. Thalmeier, 2 and T. Chatterji 3 

1 Max- Planck- Institute for the Physics of Complex Systems, Nothnitzer Str. 38, 01187 Dresden, Germany 
2 Max-Planck-Institute for the Chemical Physics of Solids, Nothnitzer Str. 40, 01187 Dresden, Germany 
3 Institute Laue-Langevm, BP 156, 38042 Grenoble Cedex 9, France 

We have studied the volume collapse of LaMnOs at the Jahn- Teller (JT) transition temperature 
Tjt = 750 K which has recently been found in high temperature powder x- ray and neutron 
diffraction experiments. We construct a model Hamiltonian involving the pseudospin of Mn 3+ e g 
states, the staggered JT distortion and the volume strain coordinate. We show that the anharmonic 
coupling between these primary and secondary order parameters leads to the first order JT phase 
transition associated with a comparatively large reduction of the unit cell volume of AV/V~ 10 -2 . 
We explain the temperature dependence of JT distortions and volume strain and discuss the volume 
change as function of the anharmonic coupling constant. A continuous change to a second order 
transition as function of model parameters is obtained. This behaviour is also observed under Ba 
doping. 

PACS numbers: 75.47.Lx, 71.70.Ej 



I. INTRODUCTION 

The parent compound of colossal magnetoresistive 
manganites, LaMnOs has drawn a lot of attention be- 
cause of its various states of spin and orbital order under 
variation of temperature 0, 0, 0] • This compound un- 
dergoes a structural phase transition at Tjt = 750 K 
associated with an orbital order-disorder transition ]4$. 
It has been observed that in the low temperature phase 
the orbital ordering is of C type with alternate arrange- 
ment of d3 X 2_ r 2 and d^ y i_ T i orbitals in the a6-plane while 
the planes are stacked along the c— axis |5j. On further 
lowering the temperature LaMn03 undergoes a magnetic 
transition to an A-type antifcrromagnetic phase at the 
temperature Tn — 145 K where the spins are aligned 
parallel to each other in the ab— plane and antiparallel 
along the c— axis 0,0]. 

In a recent experiment T. Chatterji et al. Q have in- 
vestigated in detail the Jahn- Teller transition in LaMnOs 
using high temperature x-ray and neutron diffraction on 
powder samples. They observed that the unit cell vol- 
ume of LaMnOs decreases with increasing temperature 
in a narrow temperature range below Tjt and then un- 
dergoes a sudden collapse at Tjt- It was argued that this 
striking volume collapse is caused by the orbital order- 
disorder transition. In the orbitally ordered phase the 
packing of MnOg octahedra needs more space than in 
the disordered phase. The authors compared this unusual 
phenomenon to the melting of ice which is accompanied 
by a similar volume collapse. They also observed that 
with very small doping of Ba the first order like transi- 
tion becomes second order [jj . The above mentioned ex- 
perimental observation of the volume contraction at the 
Jahn- Teller transition temperature in LaMnOs has moti- 
vated us to study the effect from a microscopic model for 
LaMnOs. In addition to the usual Jahn- Teller interac- 
tion terms we include the coupling between volume strain 
and the Jahn- Teller distortions. This allows us to explain 
both the observed volume contraction and the first order 



transitions in the temperature dependence of Jahn- Teller 
distortions Q2 and Q3 within the same context. 



II. PSEUDO SPIN MODEL FOR THE JT- 
TRANSITION IN LaMnOs 



In LaMnOs the Mn +3 ion in each Mn06 octahedra 
is in its t^gCg state, the degenerate e g orbital is singly 
occupied and hence JT-active. The local JT-distortions 
around each Mn +3 ion interact with each other coopera- 
tively and give rise to the observed orbital ordering. Sev- 
eral theoretical investigations have been reported in the 
past to study the cooperative JT-phcnomcna in 3d- com- 
pounds like MnF3, KCUF3, LaMnOs and others. Among 
the first few attempts to study such systems is the ap- 
proach of Kanamori . He started from a microscopic 
Hamiltonian with the couplings of JT- electrons to uni- 
form bulk distortions as well as to all vibrational modes. 
In this model two possible ordered phases are consid- 
ered, the 'ferroorbitaP phase where the local distortions 
together with the orbitals of JT-ions align in the same 
direction and the the 'antiferroorbitaP case where they 
align in opposite directions leading to a staggered or- 
der. Kanamori's theory was successfully applied to some 
spinel-type JT-crystals with ferroorbital order to explain 
the bulk tetragonal distortion and the nature of the struc- 
tural transition ^l|. The theory of this type of ordering 
was further developed by Pytte ^2 f° r more than one 
JT-ion in an unit cell. A review of JT phase transitions 
was given by Sturge |13| In addition it was predicted that 
the structural phase transition for the antiferro-orbitally 
ordered state should be of second order even in the pres- 
ence of higher order terms in the JT- energy and anhar- 
monic terms in the JT distortions [l(| • But the observed 
phase transitions in most of these antiferroorbitally or- 
dered compounds are in fact first order like. Later on 
Kataoka [14J showed that for large higher order JT inter- 
actions and if both ferro- and antiferro- distortions are 
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taken into account the transition may still be of first or- 
der. We will now use a similar model, supplemented by 
the effect of the volume strain, to the present case of 
LaMnC>3 . As in 0, we start from a microscopic 
model Hamiltonian with the usual first order JT-terms 
coupling the electronic states to the (E s ) type JT- 
distortions. Since the fully symmetric Ti (Ai g ) volume 
strain can always couple to the square of the JT order pa- 
rameter, irrespective whether it is of ferro- or antiferro- 
type, one has to expect a corresponding third order an- 
harmonic term in the lattice energy. It has two conse- 
quences: Firstly it induces a spontaneous volume strain 
Qo below Tjt as secondary order parameter in addition 
to the primary staggered JT order parameter Q s . Sec- 
ondly, if the coupling term is large enough compared to 
the energy of the JT distortion, it can also lead to a first 
order JT transition. Since in LaMnC>3 we indeed observe 
a large volume strain below Tjt and the transition is of 
first order, despite having a pure staggered order param- 
eter, we conclude that this term which couples volume 
and JT- strains is more important than the higher order 
terms in the JT energy itself introduced in 0, 0, ^| . 
Our model Hamiltonian is then given by 

H = \gt{Q\ + Q%) -go^2{Q2<J Xi + Q 3 a Z i) 

i 

+\c B Ql-iiQ {Ql + Ql)-ioQl (l) 

Here the first term is the elastic energy of JT dis- 
tortions with Cij — Qcij where Q is the unit cell vol- 
ume and for the ferro- distortive the bulk 
elastic constants. Thus Ct/^ = C11-C12 is the elastic 
constant of T3 symmetry distortions (Q2,Qs). In the 
ferro- distortive case they are related to the cartesian 
elastic strain components en by (Q2,Q3) = 1/V2((en 
-e22),l/v / 3(2e33-e22-en)). The second term in Eq. 
represents the first order JT-coupling to the distortions 
with go being the coupling constant. The third term gives 
the energy due to volume change where the bulk modulus 
Cb = l/3(Cu + 2C12) and Qo is the volume strain given 
by Qo — e n + e 22 + 633. Here we assume that Qo(T) 
= (V(T)-V (T))/V (Tjt) is the dimensionless volume 
change caused by the JT ordering in addition to the al- 
ready present backround volume variation of Vo (T) . The 
coupling between volume strain and the JT- distortions 
is given by the fourth term with a coupling strength 71. 
The last term represents the anharmonic energy due to 
the volume change. Furthermore <j x and a z are pseudo 
spin Pauli matrices within the subspace of e± g states. The 
experimental observations Q show that the Mn-0 bond 
lengths alternate between short and long in the xy-plane 
and are the same along z-axis which suggests that pure 
staggered orbital order is realised in LaMnC>3 without 
any admixture of ferrodistortion . The pure staggered 
order is of C-type corresponding to wave vector Q = 
(7r,7r,0). To study this case we divide the lattice into two 
sublattices A and B and the staggered order requires that 
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FIG. 1: The free energy F(Q 3 ,T) per Mn+ 3 ion is plotted 
against the staggered JT-distortion Q 3 (the primary order 
parameter) for different temperatures above and below the 
transition temperature. The value of the order parameter for 
which the system has minimum free energy is shown by solid 
dots for a sequence of temperatures starting at T/T jt = 0.6 
(deepest minimum on the right side) increasing in steps of 
0.05 to a value of T/Tjt = 1.1. At Tjt the order parameter 
Q s jumps to zero from a finite value. 



distortions in A and B have equal magnitude but oppo- 
site signs, i.e. Q 2 a = —Q2B and Q 3 a = —Qzb)- Wc 
write Q2 and Q3 in polar coordinates as Q2 = Qsin 9 and 
Q3 = Qcosfl. We then have Qa = -Qs = Qs or Q s = 

1/2(Qa-Qb). 

Thus in terms of Q s Eq. QJ then becomes, 

H = \c T Ql - goQ s ^ (sin flcr-ri + cos da Z i) 

i 

+\c B Ql - liQvQl ~ loQl (2) 

Note that for the staggered case Cr is not identical to 
the bulk elastic constant and go is different from the fer- 
roorbital case, however for simplicity we do not introduce 
new notations. The Hamiltonian in Eq. [21is invariant un- 
der rotations in the doubly degenerate e\ g subspace, the 
appropriate ground state orbitals d 3a .2_ r ,2 and c?3j,2_ r 2 on 
A- and B- sublattices corresponding to 9 = 30° at low 
temperatures are selected by higher order JT interactions 
which we do not include explicitly, for the volume con- 
traction effect this is indeed not necessary since the term 
proportional to 71 is independent of the mixing angle 9. 
The Hamiltonian is then diagonalised and the free energy 
of the system is given by 

F = -Nk B T\ogZ{Q s ) + NE 
Z(Q S ) = 2 cosh {g Qs/k B T) (3) 
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Eo = ^CtQ 2 s + ^C b Q 2 o-1iQ Q 2 s ~1oQ 3 o 

In the following we will consider the JT-distortion Q s 
as the primary order parameter of the transition and 
the volume strain, Q , as the secondary order param- 
eter. Minimising the free energy with respect to Qo we 
arrive at the relation 

Qo = XQl for l^Ql^l (4) 

where we defined A = 71 /Cb- As shown below the 
condition in Eq. (0J is well fulfilled. Substituting Eq. @ 
into the expression of the free energy and minimising it 
with respect to Q s we get the mean-field equation for the 
staggered order parameter: 

C T Qs - 2\ 2 C B Q 3 S - 6 7o A 3 Q s 5 = g tanh (g Q Qs/k B T) (5) 

For the simple case A = (no coupling of JT distortion 
and volume strain) we have a second order transition with 

k B T JT = $- F(0) = -~k B T JT 
QM = p- QoW^AQ^O) 2 ^-^ (6) 

For the general case A > the solution for the tran- 
sition temperature Tjt and the primary Q s (T) and sec- 
ondary (Eq. (@J) Qo(T) for T < T JT has to be obtained 
numerically. The evolution of the free energy and the 
associated minimum described by the solution of Eq. [S] 
is shown in Fig. ^ For the parameters chosen the first 
order nature of the transition due to the 71 (A) coupling 
term is clearly obvious. 

III. NUMERICAL RESULTS AND DISCUSSION 

We now solve Eq. for the JT distortion numerically 
to obtain the temperature dependent order parameters . 
The energies per volume are given in the unit of GPa for 
convenience and the unit cell volume is taken from to 
be V = 245. 64A 3 . Assuming that C B and Ct are known 
from measurements (this is only approximately true for 
Or) we have three fit parameters, the JT coupling go 
and the anharmonic constants 70 and 71 , or equivalently 
A= 71/CB. We fit these parameters consecutively. First 
the JT-coupling strength go = 1.898 GPa is determined 
from the experimentally observed JT transition temper- 
ature Tjt — 750 K y| using the approximate relation 
(Eq. 10) go/Cr = k B TjT valid for the second order 
transition. Then the anharmonic coupling parameters 
A (dimensionless) and 70 were determined to be 13.0 and 
-400.0 GPa respectively by fitting the theoretical curves 
of Q S (T) and Qo(T), including the jumps at Tjt, self- 
consistently to the experimentally observed temperature 
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FIG. 2: Theoretical curve for the temperature dependence 
of normalised JT-distortion (Qs/Qs(0)) (solid line) with pa- 
rameters C T = 85.0 GPa, A = 13.0, 70 = -400 GPa. The 
bulk modulus Cb(T) is taken from fl5j| (see inset of Fig- El • 
The solid dots represents the experimental data |8J for the 
temperature dependence of Q 3 . 




FIG. 3: Calculated temperature dependence of normalised 
volume strain (Qo/Qo(0.63Tjt)) (solid line) for the same set 
of parameters as in Fig.^along with the experimental data || 
(solid dots) for the temperature dependence of the normalised 
unit cell volume in LaMnOs. The temperature dependence of 
Cb is extracted from [15| (see inset). 



dependence of of JT- distortion and the unit cell volume 
□ 

Thereby we have also included the temperature depen- 
dence of the bulk modulus Cb(T). It may be extracted 
from experimental results |15| obtained for Sr- doped 
LaMnOa by properly rescaling Tjt and fitting to the 
data in the temperature range between JT- and mag- 
netic transitions. The resulting Cb(T) is shown in the 
inset of Fig. [21 In Fig. [21 we show the temperature de- 
pendence of the staggered JT-distortion Q s (normalised 
to its value at zero temperature) from our calculation 
along with the corresponding experimental points. The 
temperature dependence of the volume strain Qo (also 
normalised by the value at zero temperature) calculated 
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FIG. 4: The temperature dependence of Qo normalised by its 
value at T = is shown for two different values of A. The 
curves for A = 13 (solid line) and A = 7 (dashed line) show 
first and second order transitions respectively. The inset in 
(a) gives the dependence of the jump in Qo at the transition 
temperature on A. In (b) we show the corresponding curves 
for the temperature dependence of entropy normalised at Tjt 
for the same set of values of A. 



from our model along with the corresponding experimen- 
tal points for the unit cell volume is presented in Fig. 3. 
With the parameter set used for Figs. 12131 and using the 
relations in Eq.HJwe get an estimate -tQt-Q? — 0.9 • 10~ 2 
for the condition in Eq.|4|in agreement with the assump- 
tion. The parameters have been chosen to lead to the 
experimentally observed first order phase transition. For 
illustration we show in Fig. 4(a) that a reduction of the 
anharmonic coupling of primary (Q s ) and secondary (Q ) 
order parameters leads to a change from first to second 
order phase transition as witnessed by the vanishing jump 
in the order parameters and the associated vanishing of 
the latent heat. 

In the inset of Fig. 4(a) we show that the vol- 
ume jump actually first increases when A is reduced 
and then suddenly drops to zero. This transition from 



first to second order is indeed observed for Ba- doped 
Lai-^Ba^MnOa |j| where the volume jump vanishes al- 
ready for x = 0.025. In Fig. 4(b) we show the tempera- 
ture dependence of the entropy changing discontinuously 
(continuously) with temperature at the transition for first 
(second) order transition. 

In this work we have given an explanation to the exper- 
imentally observed first order JT- transition in LaMnC>3. 
The first order transition is due to a strong coupling of 
primary JT-distortion and secondary volume strain or- 
der parameters. This result differs from those of previ- 
ous models |lfjj | which did not include the coupling to 
the volume strain but rather consider higher order terms 
in the JT energy. There the temperature dependence of 
JT-distortions is of second order type in the pure anti- 
ferro orbital order. Therefore the earlier models cannot 
be directly applied to the present case of LaMnC>3 where 
one has a first order transition despite having a pure stag- 
gered order parameter. ^From our calculation we see that 
the term which couples volume strain and JT-distortion 
in the model Hamiltonian is also responsible for the sharp 
collapse of the bulk volume. In the absence of such cou- 
pling 0, ^| , the JT-distortions are volume conserving. 

We have not considered the spin degrees of freedom in 
our model calculation as has been discussed in [lfj and 
fl7| where it is argued that spin and orbital ordering can 
interfere with each other through the orbital dependent 
superexchange interaction. Since the magnetic transition 
temperature (Tjv — 145 K) in LaMnOa is much lower 
than the orbital ordering temperature at Tjt — 750 K, 
we neglected the effect of spin correlations on the orbital 
order. 



IV. CONCLUSION 

The pronounced volume collapse in LaMnC>3 at the 
JT-transition temperature has been studied within a JT 
pseudospin model. Such models have been used pre- 
viously to describe structural phase transitions with a 
JT distorion as order parameter in spinel- type crystals 
|l0l Il2]| . To explain the additional volume contraction as 
secondary order parameter at the orbital order-disorder 
transition we have supplemented the model with an an- 
harmonic interaction between bulk volume strain and the 
JT-distortion. We have studied the evolution from sec- 
ond to first order phase transition as function of the cou- 
pling strength and give a realistic parameter set valid for 
LaMnC>3 from comparison of the observed Tjt and the 
associated jump in the JT distortion and the volume col- 
lapse. We also have shown that under suitable change of 
coupling parameters the transition changes from first to 
second order which may simulate the observed behaviour 
under small Ba- doping of the crystals. 
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